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ABSTRACT 



This thesis provides the student of combat modelling 
with a computer program for a relatively simple combat 
model that can be used in a classroom environment for study 
and analysis. The model is an aggregated, force-on- force 
ground-combat model that uses Lanchester's aimed-fire 
equations for casualty assessment. The original version 
of the model was developed in 1979 in a previous thesis, 
and shortfalls in the original model have been overcome, 
along with the addition of several enrichments to provide 
added user flexibility, A user's manual is provided to 
facilitate user access to the model from a permanent disk 
in the V/. R, Church Computer Center. 
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I. INTRODUCTION 



As the art of combat modelling becomes more advanced, 
combat modelers are continually building more and more 
complicated models. To the beginning modeler, the ability 
to understand hov/ these models operate is difficult, if 
not impossible. As a student of combat modelling, I have 
sensed a need for a relatively simple model that could be 
easily studied in a classroom environment. It was with 
this in mind thst this project was initiated. 

At present, there seems to be no simple combat model 
available that demonstrates the basics of model building 
to the beginning student. In 1979 » Smoler (a student 
at the Naval Postgraduate School) attempted to build such 
a model for his thesis research ( Smoler ,1979) . His model 
was a deterministic, force-on-force computer model that 
used Lanchester's aimed-fire equations for casualty ass- 
essment, The general scheme of his model is shown in 
Figure 1, 

Although Smoler 's model was a bold attempt at a simp- 
listic combat model, it did have some problem areas that 
warranted investigation. Some of these problems were 
first discovered during a class project in a combat models 
class (OA 4655 ) at the Naval Postgraduate School, It is 
the purpose of this thesis to study Smoler 's model and to 
undertake an enrichment program that will make the model 
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FIGURE 1. GENERAL SCHEME 
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more realistic and flexible, while maintaining it's trans- 
parency and simplicity. 

Once enhanced, the model would then be made available 
for use as an instructional tool for combat modelling 
classes. This will be accomplished through the develop- 
ment of a user's manual to facilitate the use of the model 
by students by placing the model on a permanent disk in 
the W, R, Church Computer Center where it will be easily 
accessed by any desired user. 

The remainder of this paper will discuss, in detail, 
the problems found with Smoler's model and the methods 
used to solve these problems. It will also outline some 
new features that have been incorporated into the model 
to allow user flexibility. Finally, the User's Manual 
(Appendix A) will provide all the required information 
to enable even the novice user to utilize the model. 
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II. THE ORIGINAL MODEL 



A. GENERAL DESCRIPTION 

Smoler's original model of land combat is a deter- 
ministic model that plays combat between two homogeneous 
forces, a blue force and a red force. The blue force is 
comprised of three subunits in a static defense, with each 
subunit armed with three TOW antitank missile systems. 

The red force is composed of three subunits of three tanks 
each, attacking on pre-planned routes. The battle takes 
place on the 10 x 10 Km Fulda Box that has been developed 
and used in the STAR simulation model (Wallace, 1978). Since 
the major thrust of this paper is to alter the original 
model to a form that will be easily used for classroom 
instruction, a brief discussion of the major components 
of the original model will be presented, including problem 
areas that have motivated model changes. 

B. ATTRITION PROCESS 

The attrition process in the original model utilizes 
Lanchester "aimed-fire" equations with variable attrition 
coefficients. The Lanchester equations used are fairly 
simple and will be discussed later. However, the calcu- 
lation of the attrition coefficients is of more immediate 
interest. Smoler used the Bonder-Farrell model to calcu- 
late the coefficients, the rate at which one firer 



9 



of unit i kills unit 3 targets. These are computed 

according to: 



Aij = l/ECT,^) 



where E(T. .) is the expected time for one firer of unit i 
to kill one target of unit 3 . The E(Tj^^.) is computed 
using the Bonder-Farrell formula: 



E(Ti.) = . 



(t 



m 



ti - th + (t^ + t^)/P(KIH) + 

+ t^)/P(h|m) X (d-P(h|h))/P(KlH) 
P(hlh) - P) 



+ 



v/here 

t^ = time to acquire a target 

tj = time to fire first round following acquisition 
t^ = time to fire follov/ing a hit 
t^ = time to fire following a miss 
t^ = time of flight of a round 
P = probability of a first round hit 
P(hlh) = probability of a hit following a hit 
P(hlm) = probability of a hit following a miss 
P(KIH) = probability of a kill given a hit 
This formula holds for the conditions that the hit 
probability of any round depends only on the result of the 
previous round and no accumulated damage is considered. 
Smoler also assumed that for the TOW v;eapon system, 

P(KIH) = 1.0 and P(hlm) = P(h|h) = P, thus reducing the 
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equation to: 

E(Tij) = + tj_ + t^ + (t^ + tj)(l-P)/P 

Smoler also assumed that for tanks, P(K|H) = 1.0 and 
tj = 0, reducing the formula to: 

E(T. .) = t + t, + t (1-P)/P(hlm) 

1 j 3, X in 

utilizing these equations to calculate the the 

attrition during each time step was computed using the 
Euler-Cauchy differencing equations to approximate Lan- 
chester's force-on-force attrition differential equations, 

C. UNIT LOCATION AND MOVEMENT 

Ped and blue unit locations were handled in two dif- 
ferent 'Ways in the original model. The blue locations 
were left as user inputs, v/hile the red locations had 
been pre-determined by the model builder and could not be 
altered by the user. This allov/ed flexibility of defensive 
positions, but because of the method of determining movement 
routes for the attackers, this flexibility was limited. 

The method utilized for route determination v/as, at 
best, unrealistic. For each original red location, a 
straight v;est to east route was calculated. Each route 
was divided into 40 meter intervals, since each red unit 
was assumed to move that distance during each 10 second 
time interval. This approximated an average rate of move- 
ment of 9 m.p.h. This method of route selection, although 
easy to implement, has several significant shortcomings. 

It does not allow the attacking units to utilize terrain 
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features during movement. Also, it permits attackers to 
move over terrain that is, in real life, impassible, 

D. BATTLE TERMINATION 

In any combat model, adequate battle termination rules 
must be considered. Smoler utilized two criteria to 
terminate the battle. The first of these was annihilation 
of one of the tvio forces. This criteria is reasonable in 
an expected value model like this and was adequately 
handled. This is not the case for the second termination 
rule. 

The second rule for termination is that the distance 
betv/een red and blue forces becomes too small. The use 
of this criteria for a stopping rule is fine, provided it 
is implemented properly. In the original model, a geographic 
center was calculated for each force. It v;as these center 
of mass points that were compared to determine if the 
units v;ere too close. This method can lead to problems, 
such as those that v/ere encountered during some initial 
trial runs. In one battle, the red units, during their 
advance, had eliminated one entire blue unit and had suf- 
fered the loss of one of their own units. Since the red 
units are allowed to advance as much as I50 meters apart 
and blue units were set up so they were approximately 
1300 meters apart, there was considerable distance between 
units on the battlefield. Since the two units had been 
destroyed, the center of mass computation allov/ed the red 
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units to move completely through the blue defense, while 
both sides still had forces remaining to fight with. This 
reflects a problem with the center of mass computation as 
a method for determining distances between units. 

E. DETECTION AND FIHE ALLOCATION 

Both detection and fire allocation processes are hand- 
led well in the original model and no changes have been 
made. A detailed description of these methods is contained 
in Smoler's Thesis (Smoler, 1979) . However, a brief and 
general discussion of both processes will be included here. 

The detection phenomena is modeled in two ways. First, 
a non-firing detection can occur as a result of an observer's 
random search within his designated section of responsibility. 
A 30 ° field of vision for an observer is assumed, and the 
probability that an observer is looking in the direction 
which enables him to detect a target is computed by inte- 
grating the Limicon Function over limits that are + 15° 
from the primary direction the observer is looking. The 
Limicon Function, f(9), is the following probability den- 
sity function: 

f(e) = A + B X cos(e) 

where 

D = assigned sector width/2 
A = -B X cos(D) 

B = 1 / 2(sin(D)- D x cos(D)) 

6 = primary direction observer is looking 
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Also, A and B are chosen such that 



J f(0)de = 1 

-D 

The second method of detection is a firing detection. 
This phenomena occurs when the following happens. If a 
firer's location is within ± 15 ° of an observer's primary 
direction of observation when he is firing, he is assumed 
to be detected and is added to the observer's target list. 
This models the detection of a firer by locating a firing 
signature of a weapon. 

The fire allocation process is also modeled well by 
Smoler, Since each ’firing unit is not restricted to firing 
at one target, a fraction of each firing unit is allocated 
to each target on the firer's target list. This fraction 
is determined as a function of range to each target and of 
predetermined priorities of fire that are set so each unit's 
fire priority goes to those targets to his immediate front. 
Targets to the firer's flanks are then allocated a smaller 
percent of the available firepower. Again, the detection 
and fire allocation processes have been handled well and 
no changes to these routines v/ill be made. 



III. MODEL CHANGES 



A. GENERAL 

As the previous section outlines, there were several 
major problems discovered in Smoler's original model. In 
order to make the model more flexible and classroom useable, 
several major changes were found to be necessary. This 
chapter will provide a detailed description of these 
changes, which include conditionally-deterministic attrition, 
unit locations and movement, battle termination and move- 
ment of defenders. 

B. CONDITIONALLY-DETERMINISTIC ATTRITION 

In order to introduce a stochastic process into the 
attrition computation for this model, several options were 
available. These include the use of a Markov-process to 
determine casualties or the use of random attrition-rate 
coefficients. After considering both options, and in 
keeping with the transparent nature of this model, it v/as 
decided that stochastic attrition coefficients should be 
utilized, because the attrition-rate coefficient is a 
random quantity measuring a unit's fighting ability and 
can be realized before any given battle. The follo’wing 
procedure was used to develop a method for random attrition- 
rate coefficients. 

The attrition-rate coefficient, is a measure of 

the rate a firer of unit i attrits a target of unit j. 
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This also can be interpreted as a measure of the fighting 
ability of an element of type i. It is intuitively obvious 
that this is a variable quantity that can be affected by 
many different factors, including weather, espirit-de-corps, 
previous engagements and leadership, to mention only a few. 
In a homogeneous model, one can attempt to capture these 
random effects by developing a distribution of initial 
fighting unit capabilities, i.e., initial for each 

unit. Since no data is readily available that captures this 
phenomena, an attempt was made to fit a distribution to 
an intuitive feeling as to how this fighting ability varies 
from one unit to the next prior to a given engagement. It 
was the author’s intuitive feeling that the Aj^^'s should be 
distributed between .3 and ,8 with the majority of the units 
being rated between ,5 and ,6. 

Utilizing the above intuition, an attempt y;as m.ade to 
"fit" a distribution to these assumed fighting levels. 
Graphically, the distribution v/ould look like Figure 2. 




Figure 2. Distribution of Initial Aj^j's 
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Initially, an attempt was made to fit a Beta distribution 
to this curve. The fit was fairly accurate in the middle 
range, but was not satisfactory in the tails. The same 
was true for a Normal approximation. It was finally 
determined that the best fit was a straight quadratic fit, 
v/here a Uniform(0, 1) input variable could be used and the 
output would yield the desired random number. The quad- 
ratic fit that was used is; 

A? . = -2U^ + 2U + .3 

Once a distribution for the attrition-rate coefficients 
had been derived, the implementation into the model was 
accomplished. Since it was assumed that the fighting 
ability of each unit is a random quantity prior to a given 
battle, it v/as only necessary to obtain a realization of 
the random variable for each unit prior to the initiation 
of the battle. This realization, A?j, is determined by 
using a random Uniform(0,l) number and the above formula, 

A new attrition-rate coefficient for each unit is then 
computed during each time step using the equation: 

{ A?. (1-r/r for OSr^r 

j- j e 6 

0 for r^ 1 r 

e 

?/here 

r^ = maximum effective range of a firer's weapon 
r = current range between firer and target 
A?j = realization of unit's fighting ability 
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This equation is utilized because it varies \vith range, 
but it also is a function of A?^. , thus creating a different 
attrition-rate curve for each unit, depending on that unit's 
fighting ability prior to the battle. Graphically, as an 
example, the attrition-rate coefficient curve for an 
of .6 and a maximum effective range, r^, of 3OOO meters 
would look as in Figure 3. 




Figure 3* Attrition-rate coefficient curve for 

A? . = ,6 and r = 3OOO m. 

1 J G 

Once the above method for determining attrition-rate 
coefficients had been selected, it was coded and included 
as a user option in the program. Once the code had been 
implemented, several runs were made utilizing different 
random number seeds to compare battle results using this 
method for attrition-rate coefficient determination. These 
runs were then compared with a run using the original 
deterministic method. The results are summarized below. 

First, the model v/as run using the deterministic attrit 
ion method. For purposes of comparison, all locations. 
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movement rates, force levels and other input variables v;ere 
held constant for all runs. Combatant attrition using the 



deterministic process as a function of time is shown in 



Figure 4* 



Force 

Levels 




Figure 4, Attrition using deterministic attrition 



In this battle, termination occurred as a result of forces 

being too close together, with force levels at termination 

reflected in Figure 4» This was then compared to the runs 

that utilized the new attrition module. The results of 

four of these runs are in Figure 5 . The A? . ' s are the 

^ 3 

realizations of the random variable denoting a unit's 
initial fighting capability prior to the battle. The three 
numbers under the red and blue headings are the realizations 
for each subunit in the battle on which casuality assessments 
are computed. 

These results show, in fact, that a randomly selected 
A? . does have an effect on the final outcome of the battle. 
This module now provides the user the option of selecting 
a conditionally-deterministic process that still uses 
Lanchester's equations as the basis for the attrition 
computation. 
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A? 


Force Level 




Blue 


Red 


Blue 


Red 


Time 


Result 


.33 .77 .46 


.79 .49 .56 


1.7 


1.5 


780 


Too close 


.39 .67 .58 


.63 .78 ,36 


4.3 


0 


330 


Blue vd.ns 


.78 ,67 .80 


,74 .68 .72 


3.2 


0 


320 


Blue vn.ns 


.32 ,43 .68 


,51 .77 .65 


1,2 


3.8 


640 


Too close 


Deterministic 


1,1 


2.8 


460 


Too close 



Figure 5« Comparison of Stochastic and Deterministic 
Battles 

C. UNIT LOCATION AND MOVEMENT 

As was pointed out in an earlier section, the unit 
locations and movement logic that Smoler used tended to 
be rigid and inflexible. The changes that have been made 
to rectify this and allow more flexibility for the user 
will be discussed here. 

First, the user is now responsible for inputing all 
combatant locations, both attackers and defenders. The 
format for these inputs will be specified in Appendix A, 
the User’s Manual, Aditionally, the user is given the 
option of selecting routes for the red advance or allowing 
the model to calculate straight west to east routes as in 
the past. To add more flexibility, the user vn.ll also be 
required to input the rate of advance (vehicle speed) for 
the attacking forces. The inputed rate of advance is 
used in both methods of route calculation, thus controlling 
the* red attack. 
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1, straight Line Routes 



The option to use straight line routes from west 
to east is handled in the same manner as before. The route 
is broken down into discrete distance intervals from the 
initial red location straight to the east. The only dif- 
ference is that the length of each interval is determined 
by the vehicle speed that the user inputs. The route for 
each red unit consists of 125 equal length intervals. 

2. User Determined Routes 

The option to allow the user to select routes for 
the red advance has been added to enable the user to make 
use of available terrain and to add realism to the model. 

The method for calculating these routes is straightforward 
and is discussed next. 

The user is required to input the original location 
of each red subunit and from one to ten nodes he wishes 
each attacking subunit to move through. This information, 
along y/ith vehicle speed, is used to calculate route inter- 
vals that move the attacking unit through each of the 
designated nodes, A complete route would look like the one 
depicted in Figure 6, The method used to compute the routes 
is as follov/s. 

The straight line ground distance between the first 
two adjacent nodes, DIST, is calculated as shown in Figure 
7, The angle between the desired direction of movement 
and straight west to east movement, a, is then calculated. 
Utilizing these quantities and the distance desired to 
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Figure 6. User Determined Routes 



move during each time step, DST, the distance to be moved 
in the x and y direction, XLM and YLN, is now computed as 
shown in Figure 7. These distances are then added to the 
coordinates of the previous interval endpoint, point C 
in the figure, to determine the coordinates of the next 
interval endpoint, point D. This same distance is again 
added to compute the coordinates of the next endpoint, 
point S, This process is continued until the distance 
from the last endpoint computed to the next node is less 
than DST. This process is repeated between the next tv/o 
nodes until the entire route is completed. In order to 
insure all intervals are of equal length, the computation 
of the first interval between any two nodes must be con- 
sidered separately, by taking into account the distance 
left over from the last computation between the previous 
two nodes. To do this, the first interval takes the rem- 
aining distance, e, and adds it to an interval length of 
DST-e for the first interval between any tv;o nodes. This 
insures that each interval along the route is of length DST, 
which is the required length. 

D. BATTLE TERMINATION 

As previously stated, the battle is term.inated either 
by annihilation of one of the two forces or by forces 
getting too close. The annihilation criteria has not been 
changed. However, to insure the distance criteria is 
effective, two changes have been made. 

First, the distance between units is no longer calculated 
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NODE COORDINATES 

A XIC(N-1 ) ,YIC(N-1 ) 

B XLOC(I,J) 

C XIC(N) ,YIC(N) 

D XIC(N+1 ),YIC(N+1 ) 

E XIC(N+2) ,YIC(N+2) 

F XLOC(I,J+1) 



USE 

(N-1)^^ Interval Endpoint 

User Inputed Node 

N^^ Interval Endpoint 

Interval Endpoint 

(N+2)’^'^ Interval Endpoint 
User Inputed Node 



DIST = NXL^ + YL^ 
a = Tan"^ (Y/XL), where Y = 1yl| 

e = distance less than DST at end of calculation 
of intervals between adjacent nodes. 

YIII = DST X SIN(a) XLN = DST x COS(a) 



Figure 7. Route Computation 






as a center of mass distance. Instead, a distance is 
computed between each attacking subunit on v/hich casualties 
are being computed that is still alive and each defending 
subunit that is still in the battle. If any of these 
distances becomes too small, the battle is considered to 
have moved to a "close-in" combat (hand-to-hand) that this 
program does not currently model. For this reason, the 
battle is then terminated. However, to insure the attackers 
do not pass by the flanks of remaining defenders and remain 
outside termination distance, a check of x coordinates 
for each subunit is also made. If any attacking subunit's 
X coordinate places him past the location of the forward 
most blue defensive subunit still in the battle, the battle 
is also terminated. 

The criteria for being too close has been left as a 
user determination. It is one of the inputs that is 
described in the User's Manual, Appendix A, This allov/s 
for flexibility of breakpoint distances which also lends 
itself to the study of optimum breaking distances for 
various weapon system.s on the battlefield, 

E. ALTERNATE DEFENSIVE POSITIONS 

The last addition to the model, added to increase the 
flexibility and realism of the battle, is the option of 
having the defending units move to alternate positions if 
the attackers close to within a user specified range. 

The move is handled in a simple and transparent manner. 

If the user selects this option, the defender will move 



25 



v/hen the attacking forces close within the breakpoint 
distance. V/hen this occurs, each defensive unit that is 
still alive v/ill move to an alternate position the user 
has selected. The duration of the move is also a user 
input. He simply specifies the number of 10 second time 
steps he wants to allow the blue forces to move. At the 
completion of the move, during which the red forces continue 
to attack, the battle will continue until one of the term- 
ination criteria discussed above has been met. 
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IV. FUTURE ENHANCEMENT 



Although the above changes have been implemented and 
a more useable and flexible model has been created, there 
are several areas that could still be considered in future 
v/ork. These include the introduction of heterogeneous 
forces, explicit computation of ammunition expenditure 
and artillery (indirect fire). Possible methods for 
employing these ideas are discussed below, 

A. HETEROGEJrEOUS FORCES 

- - — The current model involves combat between homogeneous 
forces only. In other words, each opposing force is 
comprised entirely of one weapon system type. Added 
flexibility could be attained by allowing heterogeneous 
force structures on both sides. This would enable the user 
to investigate the effect of different force mixes on 
battle outcome. 

The introduction of different weapon systems within 
a single unit would require extensive changes to the att- 
rition process. Although Lanchester equations for aimed- 
fire could still be utilized, casualty assesment against 
each weapon system type by each opposing weapon system 
type would have to be calculated. The total attrition of 
any particular unit would then be the sum of the damage 
assessed to each weapon system of that unit, Aditionally, 
the force level for each 'weapon system would be required 
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as an input as v/ell as a separate set of hit and kill 
probabilities for each weapon system type against each 
type of target. 

Introduction of this feature would create a more comp- 
licated model, adding realism but detracting from the current 
simple and transparent form. Since the purpose of the 
current effort has been to maintain this simplicity, this 
option has not been included, but could be considered for 
future work. 

B. AiMMUNITION EXPENDITURE 

The combat process is a complicated and intricate 
process. One of the hardest, and not many times attempted, 
areas to model is the logistic area. Hov/ever, in the 
model currently being studied, the problem of ammunition 
expenditure could be modeled. 

Since the model is an expected value type model, it 
would not be effective to model the amount of ammunition 
expended by simply counting bullets as they are fired. 
Actually, in this type of model it would not be possible 
to count each round, because of the aggregated nature of 
each unit. The only way to model ammunition useage is to 
model it in the same way as attrition. This simply means 
that each unit would have to be given a starting level of 
ammunition on hand, a basic load, and each time the unit 
fired, the expected ammount of ammunition expended would 
be subtracted from what is on hand. The ammount expended 
would necessarily be a function of the size of the unit 
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firing, the number of targets fired at and the rate this 
unit is firing. This can be directly correlated to the 
attrition-rate coefficient for the firing unit. 7/hen 
ammunition on hand reaches the zero level, the unit v;ould 
have to be removed from the battle permanently or for some 
specified time period, to simulate resupply. 

C. ARTILLERY 

In any armed conflict, there are more ways to inflict 
casualties than just direct fire. Many types of indirect 
fire are utilized on the battlefield, including artillery, 
close air support and naval gunfire. To include this type 
of play in this model, primarily artillery, would involve 
a major effort with large program ’additions. 

To model artillery in an expected value model that 
already uses Lanchester equations for direct fire would 
suggest the use of Lanchester' s area-fire equations. These 
equations are no more complicated to handle than the aimed- 
fire equations that are already in use. However, other 
model considerations would have to be investigated. These 
would include location of artillery units, v/hether to use 
forward observers for target location or only use pre- 
planned strikes or both and fire allocation procedures 
to be used. Also, a significant ammount of data would need 
to be collected concerning weapon types to be used, eff- 
ective ranges, killing radii and accuracy data. 
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To model artillery ’.vould be a significant addition to 
the model. Ho^vever, this too would effect the simplicity 
and transparency the model currently possesses. This 
addition has also been left for future consideration. 
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APPENDIX A 



User’s Manual 



I. INTPODUCTION 

Smoler's model of land combat is a force-on-force 
combat model that utilizes Lanchester's aimed fire equations 
for casualty assessment. The battle is simulated on a 
10 X 10 KM piece of terrain representing an area east of 
Fulda, West Germany, A portion of the terrain map is 
inclosed as Figure 1, It is a computer model that is 
coded in FORTRAN and is available for use in the W, R. 

Church Computer Center, 

The purpose of this manual is to familiarize the user 
with the model and to provide the required Job Control 
Language (JCL) and user inputs to run the model, A 
sample listing of program output will also be provided to 
familiarize the user with the expected program output. 
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Figure 1 • Model Terrain Map 



II. PRCGRAI-1 STRUCTURE 



Smoler's enhanced model is a computerized model, coded 
in FORTRAN, containing a main program and nine subroutines. 
To assist the user in understanding hov/ the model operates, 
a brief description of the function of each of these major 
parts of the program will be included* 

A. MAIN PROGRAM 

The main program has several important functions. 

All of the input and output functions, except the line-of- 
sight data, are contained in the main program. Aditionally, 
the main program is used to structure all of the other 
functions during each 10 second time step. Attrition, 
detection, movement and fire allocation are also handled 
in the main program. The nine subroutines provide needed 
input numbers for the above calculations. The general flow 
of the program is depicted in the flow chart shown at 
Figure 2. 

B. SUBROUTINES 

There are nine subroutines included in the model. Each 
of these subroutines performs a distinct function, each of 
which will be discussed below. 

1 . Subroutin e SETUP 

This subroutine is used to read in the terrain 
data for the Fulda Map. This terrain data will be used 
when computing line-of-sight between targets and observers. 



33 



STAPT 




FIGURE 2. FLOW CHART 
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as v/ell as providing a grid system for unit locations and 
movement . 

2. Subroutine LOS 

This module, developed by Professor Jim Hartman 
of the Naval Postgraduate School, computes a percent of a 
target visible to a particular observer, given the coor- 
dinates (location) of both. This visible fraction is used 
in the detection and attrition modules in the model. 

3. Subroutine KOVEP 

This subroutine is used by subroutine LOS in deter- 
mining what portion of a particular target is covered by 
the terrain between the target and observer. This number 
is used in the detection and attrition modules. 

4» Subroutine ETK 

This module computes the expected time for a given 
firer to kill a given target. The calculation is a function 
of range, time of flight for a round and hit and kill 
probabilities for the firing weapon system. It is a number 
that is used in the comiputation of the determ.inistic attrition 
coefficients, A. .. 

-L J 

5* Subroutine STOCK 

This is used when a user has selected a stochastic 
attrition option to compute the attrition coefficients 
during each time step. The calculation is a function of 
the original stochastically determined A. . as well as a 

-L J 

function of range. 
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7. Subroutine LAMBDA 



This subroutine is used in conjunction with the 
LOS routine to compute the detection rate of a target by 
an observer given the percent of target visible to the 
observer. 

8. Subroutine ROUTE 

This is used to compute the route of each attacking 
red unit when the user has selected the option of inputing 
attacker routes. It calculates the coordinates of each 
interval endpoint along the route, making each interval 
length (distance moved during a 10 second time step) the 
same. The interval length is determined by the speed the 
user has selected and input ed for the current battle. 

9. Subroutine SLEV 

A subroutine that is used to calculate the terrain 
elevation for a given set of X, Y coordinates. This function 
is used in conjunction with the LOS subroutine in computing 
line-of-sight between observer and target. 
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AVAILABI.^ OPTIONS 



The enhanced version of this model has been written 
to allow for maximum user flexibility while maintaining 
the simple and transparent nature of the model. To allov; 
this flexibility, there have been several user options 
incorporated into the model, Each of these options, 
including user responsibilities, will be discussed here, 
with the required input data for each being outlined in 
the next section, 

A. STOCHASTIC VS DETERMINISTIC ATTRITION 

The user is required to specify whether he wants 
to use stochastic or deterministic attrition calculations. 
Both m.ethods utilize Lanchester's aimed-fire equations, 
the difference being the method of calculating the attrition- 
rate coefficients. The deterministic procedure uses the 
Eonder-Farrell method of calculating the v/hile the 

stochastic method selects an initial random for each 
unit, and uses this as a function of range to calculate 
the A^j for each time step. 

B. ATTACK ROUTES 

The second major option available to the user is 
the method of route computation for the attacking forces. 

The user has to decide whether to allow the program to 
compute straight west to east routes or to input the routes 
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he desires each attacking unit to follov;. To select his 
own routes, a user must input the number of nodes he 
wishes to have on each route and the coordinates of each 
of these nodes. The program will then compute routes 
through each node. The nodes must be inputed in order 
from west to east and should not create an angle between 
the west to east axis and the route direction that exceeds 

) cO 
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C, ALTERNATE DEFENSIVE POSITIONS 

The third option the user must consider is the 
use of alternative defensive locations. This option permits 
the user, if he desires, the capability to move the defen- 
ders to alternate positions if the attackers close ’within 
some specified distance. This breakpoint distance is 
decided and inputed by the user and is also used as the 
distance for battle termination. The alternative to moving 
the defenders is to terminate the battle when the breakpoint 
distance is reached. 

D. OTHER INPUTS 

There are other inputs that are required by the 
program that the user must provide. These include force 
sizes, weapon characteristics, unit locations and hit and 
kill probabilities. The required formats for all inputs 
are outlined in the next section. 
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IV. INPUT DATA 



All user input data is read in from cards at the start 
of the main program, MAIN. A hrief description of each 
data element as v/ell as the required format for inputing 
this data is discussed next, A sample data deck has 
also been included as Figiare 3> with the referenced data 
element number appended on the right. This sample data 
deck reflects the use of all data elements. However, 
not all the elements described below are required, and 
those that are needed only if a particular option is 
selected are noted with an asteric next to ohe data element 
number. 

Data Element Descrintion and Format 

1. Selection of attrition option and random 

number seed, A 1 is entered for deter- 
ministic or a 0 for stochastic attrition. 
This number is followed by a five digit 
random number seed. Format for this card 
is II, IX, 15. 

2, Number of blue units (NBU) and number 

of red units (NRU), Both of these 
elements are two digit integers. Card 
format is I2,iX,I2. 

3* Effective v/eapon ranges listed in the 

following order. Red minimum, red 
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4 . 



5. 



6 . 



7. 



* 



maximum, blue minimum and blue maximum. 
Card format is if(F6.1,lX) . 

Red force levels for each unit. 

Card format is F3»l) IX. The card v/ill 
contain the number of entries that equals 
the number of red units (^IRU) in data 
element 2. 

Type of route com.putation desired followed 
by vehicle speed. Use a 1 for user 
determined routes or a 0 for program 
determined routes. A one digit entry will 
designate desired speed as follov/s: 

1 for 9 ‘Tiph 

2 for 12 mph 

3 for 15 mph 

4 for 13 mph 

Card format for these elements is 12, 

IX, 12. 

X, Y coordinates for each initial red 
location. One card is needed for each 
red unit with format F6.1, IX, F6.1. 

If the user has selected to enter his 
own routes in data element the 
following is required for each route, 

a. The number of nodes in the route 
(from 01 to 10), This card's format 
is 12. 
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b. X, Y coordinates of each node 
along the route in order from west 
to east. Cne card for each node in 
format F6.1, IX, F6.1. 

X, Y coordinates for blue location, force 
level for each unit, primary search direction 
for that unit and desired search width. 

Force level must be betv;een 1.0 and 3»0j 
search direction betv/een 135 ° and 225 °, 
and search width from 30° to 120*^. One 
card for each blue unit with format F6.1, 

IX, F6.1, IX, F3.1, IX, 13, IX, 13. 

Specify if blue is to move to alternate 
defensive positions, breakpoint distance, 
and num»ber of 10 second time intervals 
allowed for the move. A 1 for no move 
or 0 for option to move. Format for this 
card is II, IX, F6.1, IX, 12. 

Alternate blue position X and Y coordinates 
if move is desired (data element 9). One 
card for each location v/ith format F6.1, 

IX, F6.1 

Red weapon system hit and kill probab- 
ilities. The probabilities are entered 
for 6 range bands, with one card for each 

range band. Four probabilities for each 

s t 

range are probability of a 1 round hit. 



probability of a hit given a previous round 
hit, probability of a hit given a previous 
round miss, and probability of a kill 
given a hit. The 6 range bands are from 
0 to 500, 500 to 1000, 1000 to 1500, 

1500 to 2000, 2000 to 2500 and 2500 to 3OOO. 
Each of these six cards has the format 
4(F4.2, IX). 

12. Same as item 11, except the probabilities 

for the blue v/eapon systems are entered. 

Six cards with format 4(F4.2, IX), 

In the formats listed above, a IX is a space, an Ix is an 
integer of x digits, and a Fx.y is a real number of length 
X with y digits to the right of the decimal. For example, 
an 12 could be a 25> an F6.1 could be 3487.4* A format 
containing these specifications could be 12, IX, F6.1. 

The data card would then look like 25 5487. 4> v/ith the 25 
beginning in card column one. 



card 1234 ifiMent 

column 12345678901234 % 789012345678901234%7890 number 



1 28945 


1 




03 03 


2 




0000.0 2500.0 0500.0 4000.0 


3 




3.0 3.0 3.0 


4 




1 2 


5 




2500.0 7300.0 

2400.0 5900.0 

2450.0 5300.0 


6 




02 

4900.0 7200.0 

7100.0 6700.0 
02 

5350.0 6000.0 

7000.0 6300.0 
01 

6900.0 5500.0 


7 




5700.0 7800.0 3.0 190 120 

5800.0 6500.0 3.0 180 120 

5900.0 5900.0 3.0 180 120 


8 




0 0500.0 4 


9 




7000.0 7200.0 

7100.0 6500.0 
7100.0 5800.0 


10 




0.85 0.85 0.75 0.70 
0.80 0.80 0.75 0.70 
0.75 0.75 0.70 0.65 


11 





0.60 0.65 0.60 0.55 

0.45 0.50 0.50 0.35 
0.20 0.20 0.20 0.20 



Figure 3 « Sample Data Deck 
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card 

column 



1 3 4 

1234567890123436789012345678901234567890 



data 

element 

number 



0.60 0.70 0.65 0.85 
0.85 0.90 0.85 0.90 
0.80 0.85 0.85 0.80 
0.75 0.80 0.75 0.70 
0.60 0.70 0.65 0.65 
0.40 0.45 0.40 0.50 



Figure 3» Sample Data Deck (cent) 
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V. REQUIRED JCL 



The model and line-of-sight data are residing on 
permanent disk in the W. R, Church Computer Center. In 
order to exercise the model, the required job control 
language (JCL) is illustrated in Figure 4* Due to the 
CPU time required for execution, a time parameter of 
TIMS = 2 must be used on the job card. 

//JOB CARD , TIME=2 
// EXEC FORTCLG 

//FORT.SYSIN DD UNIT=2314, V0L=SER=PAT002,DISP=SHR , 

// DSN=S1360.SM0LER . PGM , DCB= (RECFM=FB, LRECL=80 , BLICSIZE=3200) 
/* 

//G0.FT08F001 DD UMIT=2314,V0I^SEE=PAT002,DISP=SHR , 

// DSn=S1360.L0SDATA,DCB=(RECFM=FB, LREC I^80,BLKSIZE=3200 ) 

/^ 

//GO.SYSIN DD 

USER DATA DECK 



/^ 



Figure Required JCL 
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VI . EXPECTED OUTPUT 



Once the model has been executed, the user can expect 
the following output: 

1 - Program listing of the model. 

2 - A summary of the initial battle conditions, 

including starting locations and options selected, 

3 - A summary of battle conditions after each 10 

second time step, including unit locations, force 
levels, unit status, percentage lost and targets 
on each unit's target list. The status will 
show one of the following: 

0 - Unit alive, not firing 

1 - Unit alive and firing 

2 - Unit killed 

3 - Unit moving 

if - A statement of the reason for battle termination, 

A sample of the initial battle summary is at Figure 5 and 
the time step summary is at Figure 6. 
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INITIAL BATTLE INFORMATION 



UNIT 


X 


Y 


FORCE LEVEL 


1 


1900.0 


7800.0 


6.0 


2 


1700.0 


6500.0 


6.0 


3 


1700.0 


4900.0 


6.0 


4 


5800.0 


7700.0 


3.0 


3 


7400.0 


6300.0 


3.0 


6 


6200.0 


5150.0 


3.0 



ATTRITION IS STOCHASTIC 



ROUTES DETERMINED BY USER 

RED VEHICLE SPEED IS 15.0 

BREAKPOINT DISTANCE IS 1500. 0 

BLUE WILL MOVE TO ALTERNATE POSITIONS 
ALTERNATE POSITIONS ARE: 

UNIT X Y 

k 9300.0 7300.0 

5 9800.0 6200.0 

6 9200.0 5800.0 



RED KILL PROBABILITIES 



RANGE 


p 


PHH 


PHM 


PKH 


500 


0.95 


0.95 


0.95 


1 .00 


1000 


0.80 


0.83 


0.81 


0.95 


1500 


0.65 


0.70 


0.70 


0.90 


2000 


0.40 


0.50 


0.55 


0.75 


2500 


0.30 


0.40 


0.40 


0.75 


3000 


0.20 


0.30 


0.25 


0.60 


BLUE KILL 


PROBABILITIES 


RANGE 


P 


PHH 


PHM 


PKH 


500 


0.80 


0.80 


0.80 


1 .00 


1000 


0.80 


0.80 


0.80 


0.95 


1500 


0.70 


0.70 


0.70 


0.90 


2000 


0.65 


0.65 


0.65 


0.85 


2500 


0.45 


0.55 


0.45 


0.70 


3000 


0.45 


0.55 


0.45 


0.60 



Figure 5. Initial Battle Summary 
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TIME = if 80 SECONDS 



UNIT 


X 


Y 


FORCE LEVEL 


STATUS 


LOST-PCT 


TARGETS 


1 


4933.2 


7193.2 


5.0 


0 


0,160 




2 


4832.3 


5799.9 


4.7 


0 


0.214 




3 


4902.6 


4579.8 


6.0 


0 


0.0 




4 


5800.0 


7700.0 


0.0 


2 


1 .000 




5 


7400.0 


6300.0 


3.0 


1 


0.0 


2 


6 


6200,0 


5150.0 


3.0 


3 


0.0 





Figure 6. Time Step Summary 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

18 

lU 

IS 

16 

17 

18 

IS 

20 

21 

22 

29 

24 

25 

26 

27 

28 

29 

30 

91 

32 

99 

94 

95 

96 

37 

38 

39 

40 

41 

42 

49 

44 

45 

46 

47 

48 

49 

SO 
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APPENDIX 3 

Program Listing 



COHHON /GRPl/ IPRDIR (6) . ISECHD (6) .HVTOIR (6) ,X (6) .T (6) .SPQ (6) 
COMMON /GRP2/ Tfl (2J ,T1 (2) ,TH(2) ,TMt2l .TFl (2) ,TF2(2) ,TF3 (2) , 

IP 12.6) ,PMH (2.6) .PHM (2.8) .PKH (2.6) .TF (2) 

COMMON /GRP3/ NBU.NRU.FL (6) .FO (6) .NOl (3) .XIC (3.200) .TIC (3.200) . 
1I01R0.200) .flVSP.lSPO 

I.IUSTAT (6) . 11 (6) .LOST (6.6) . VlSFRfl. VISFRB.SIZETK, 

131ZETW.NT (6) .NF (6) .SRF.DlSMflX. 

INLOSC (6.6) .VISFR (6.6) .RMINTK.RMXTK.RMINTH.RMXTW.OP.TOHFR.TNKFR. 
IPTT (3.3) .RF.POfl(6.6) .flP0fl(6.6) .LOR (6.6) .NR (6) .OFL (6) .POL (61 
COMMON /GRPU/ TPOL (6) .OLDQ (6,6) .0 (6.6) 

COMMON /GRP5/ LOT (6. 6) . ROT (6. 6) 

COMMON /HILLS/ XC (100) . TC (100) .PERK (100) .SX (100) .ST (1001 .RHOClOOl 
COMMON /HILLS/ SCALE ( 1 00) . TMORHO (1 00) . TUOSCL ( 100) . BASE 
COHHON /HILLS/ NHlLLS 

COMMON /COVER/ CXC (150) .CTC (150) .CPEAK (150) .CPXX (ISO) .CPTT (150) 

COMMON /COVER/ CPXT (150) .NCVELS 

COMMON /COUNTR/ KH.KHH. KV.KN.KGRS.KELL.KINT 

COMMON /GRID/ LST (10. 10) .NHL (10. 10) .LISTH (M50) .KHREP (100) .KTREP 
COMMON /GRID/ LSTC (10. 10) .NC (10. 10) .LISTC mOOl .KCREP (150) 

COMMON /GRP6/ ALPHA (6) 

COMMON /GRP7/ XR (6) . TR (6) . IHOVE (6) 

INITIRLIZRTION. 

CALL OVFLOH 
8L»0. 0 
RL-0.0 
MP»0 

Pfll-22. 0/7.0 
ZL«.00Q01 
CftLL SETUP 

REflO TERRAIN DATA FOR LINE OF SIGHT 
CHECK FOR STOCHASTIC OR DETERHINISTIC ATTRITION 

REA0(5,130) ITRIT.IS 
130 FORMATdl.lX.lS) 

00 192 I«K6 

CALL RANDOM (IS, TRAN. n 

ALPHA (1) « (-2.«YRAN«m2) + (2. kTRAN*^ . 3) 

132 CONTINUE 



READ IN NUMBER OF BLUE AND RED UNITS 

READ (5,200) NBU.NRU 
200 FORMAT (12. IX, 12) 

INITIALIZE HEAPQN SIZES 

SIZETK-2.5 



k9 



51 

52 

53 

5M 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 



o o 



SlZETW-2.5 



C 

C READ IN EFFECTIVE WEAPON RANGES 
C 

READ (5,102) RMINTK,RMXTK,RMINTW,RMXTW 

102 FORMAT (F6. 1, IX, F6. 1 , 1 X, F6. 1 , 1 X, F6. 1 , 1 X) 
C 

C INITIALIZE PM,RF,TOWFR,TNKFR AND NOD 
C 

PM=.352 
RF». 5 
T0WFR».03 
TNKFR*. 1 
NOD-2 

DO 101 I-KNRU 
NOI (1) =125 
101 CONTINUE 
K*NRU+1 
L-NRU-^NBU 
DO 111 1»1,L • 

11 ( 1 ) -0 
111 CONTINUE 
C 

C READ IN FORCE LEVELS OF EACH RED UNIT 
C 

READ (5, 103) (FL (I) , I»1 ,NRU) 

103 FORMAT (3 (F3. 1 , IX) ) 

CHECK FOR TYPE OF ROUTE DETERMI N I T I ON 

READ (5, 106) IRTE, ISPD 

106 F0RMAT(I1,1X,11) 

IF (ISPD.EQ. 1) AVSP=9.0 
IF (ISPD.EQ. 1) DST-40.232 
IF (ISPD.EQ.2) AVSP-12.0 
IF (ISPD.EQ. 2) DST-53.643 
IF (ISPD.EQ. 3) AVSP-15.0 
IF (ISPD.EQ. 3) DST=67.053 
IF (ISPD.EQ. 4) AVSP-ia.O 
IF (ISPD.EQ. 4) DST=80.463 

C READ IN INITIAL RED LOCATIONS 
C 

DO 6 1=1, NRU 

READ (5, 107) XIC(I,1),TIC(I,1) 

107 F0RMAT(F6.1,1X,F6.1) 

6 CONTINUE 

IF (IRTE.EQ. 1) 60 TO 108 
DO 2 1=1, NRU 
DO 2 J=2,125 
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II » 



101 

102 

103 

loy 

105 

106 

107 

108 

109 

no 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 



TIC (I, J) »Y1C (I . J-1) +DSTh (J-1) 
xicn, j) *xic (I, J- 1 ) •►dstm (j-n 
lOIRd, J) »0 

2 CONTINUE 
GO TO 109 

108 CALL ROUTE 

109 SUMR0»0.0 

00 3 I»1,NRU 
FOdl-Fim 
SUMRO-SUMRO'^FO (I) 

X d) =XIC (1, 1) 

Yd)=Yicd.n 

MVTDIRd) »IDIRd, n 
SPDd) •flVSP 
lUSTflT m »0 
IPROIRd)*IDIRdn) 

ISECHDd) »120 
NF d) -1 
II d) >1 

3 CONTINUE 

c 

c READ IN INITIAL BLUE LOCATIONS 
C 

SUMBO-0.0 
DO 4 I»K,L 

READ (5,104) X (I) , Y d) ,FL (I) , IPRDIR (I) . ISECWO (I) 

104 FORMAT (F6. 1 , IX, F6. 1 . 1 X, F3. 1 , IX, 1 3, 1 X, 1 3) 

FO d) =FL d) 

SUMB0»SUMB0'^F0 (I) 

MVTDIR d) «0 
SPO (I) *0.0 
lUSTATd) *0 
IMOVE d) *0 

4 CONTINUE 
C 

C CHECK FOR ALTERNATE BLUE POSITIONS AND READ IN IF WANTED 
C 

READ (5,400) lALT, BREAK, ITEM 

400 F0RMAT(I1,1X,F6.1,1X,I2) 

IF (lALT.EQ. 1) GO TO 401 
DO 402 I»K,L 

READ (5, 107) XA (I) ,YA (I) 

402 CONTINUE 

401 DELT*10. 

TA(1)=20. 

Tl (1) *8. 

TH (1) *8. 

TM (1) *10. 

TFl (1) *1. 

TF2(1)*1. 



51 



151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 
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177 

178 
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C 

C 

C 



TF3 m •!. 

Tfl(2)«20. 

T1 (2) -8. 

TH(2)«8. 

TM (2) «15. 

TFl (2) -10. 

TF2(2)»12. 

TF3(2)«15. 

READ IN HIT AND KILL PROBABILITIES 

00 5 1-1.2 
DO 514 J-1.6 

READ (5.515) P ( 1 , J) . PHH ( I , J) . PHM (I . J) . PKH (I . J) 

515 FORMAT (4 (F4.2. IX) ) 

514 CONTINUE 
5 CONTINUE 
PTT (1, 1) «1.0 
PTT (1.2) *0.8 
PTT (2.2) *0.2 
PTT (1.3) *0.8 
PTT (2,3) -0.15 _ 

PTT(3.3)-0.05‘ 

DO 31 I-l.NRU 
DO 31 J-K.L 
NLOSC (I, J)-0 
NUOSC (J. I) *0 
Q(I.J)«1.0 
Q (J. I) -1.0 
VISFRd. J)*0.0 
VISFR (J. I) -0.0 
31 CONTINUE 
IC-1 
C 

C PRINT INITIAL BATTLE INFORMATION 
C 

HRITE (6.599) 

599 FORMAT (M MX. MNITIAL BATTLE INFORMATION*) 

HRITE (6,600) 

600 FORMAT (/IX, 'UNIT*, 7X. *X*.8X, *T*.4X, 'FORCE LEVEL*) 
DO 601 I-l.L 

HRITE (6,602) I , X (I ) . Y ( 1 ) . FL ( I ) 

602 FORMAT (IX, I3,3X,F7. 1 .2X.F7. 1 . 7X.F3. 1) 

601 CONTINUE 

IF (ITRIT.EQ.l) GO TO 603 
HRITE (6,604) 

604 FORMAT (/IX, 'ATTRITION IS STOCHASTIC*/) 

GO TO 605 

603 HRITE (6,606) 

606 FORMAT (/IX, 'ATTRITION IS DETERMINISTIC*/) 
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211 
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229 

230 

231 

232 
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245 

246 

247 
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249 
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605 IF (IRTE.EQ.O) GO TO 607 
MRITE (6.608) 

608 FORMAT (IX. 'ROUTES DETERMINED BT USER'/) 

607 WRITE (6.609) RVSP 

609 FORMAT (IX. 'RED VEHICLE SPEED IS \F4.1/) 

WRITE (6.610) BREAK 

610 FORMAT (IX. 'BREAKPOINT DISTANCE IS '.F6.1/) 

IF (lALT.EQ.O) GO TO 615 

WRITE (6.620) 

620 FORMAT (IX. 'BLUE WILL NOT MOVE TO ALTERNATE POSITIONS V) 
GO TO 625 
615 WRITE (6.630) 

630 FORMAT (IX. 'BLUE WILL MOVE TO ALTERNATE PQS I T I QNS ' /IX. 
('ALTERNATE POSITIONS ARE: '/IX, 'UNIT ' .5X. 'X ' . 8X. 'T') 

DO 635 I»K.L 

WRITE (6.640) I.XA (I) ,TA (I) 

640 FORMAT (1 X. 1 3. 3X, F7. 1 , 2X . F7. 1 ) 

635 CONTINUE 
625 IRAN-500 

WRITE (6.645) 

645 FORMAT (/4X. 'RED KILL PROBABIL I T lES '/ 1 X . ' RANGE ' . 4X. ' P 
14X. 'PHH'.3X. 'PHM'.3X. 'PKH') 

DO 650 1-1.6 

WRITE (6.655) I RAN. P (1 . I ) . PHH (1 . I ) . PHM (1 . I ) . PKH (1 . 1) 

655 FORMAT (2X, 14.4 (2X.F4. 2) ) 

IRAN-IRAN^SOO 
650 CONTINUE 
IRAN-500 
WRITE (6.660) 

660 FORMAT (/4X. 'BLUE KILL PROBABILITIES '/IX. 'RANGE '.4X, 'P'. 
14X. 'PHH'.3X. 'PHM'.aX, 'PKH') 

DO 665 1-1.6 

WRITE (6.655) IRAN. P (2. I) .PHH (2.1) .PHM (2. I) .PKH (2, II 
IRAN-IRAN-^500 
665 CONTINUE 

WRITE (6.670) 

670 FORMAT ('1 '.lOX, 'BATTLE BEGINS'//) 

C 

C UPDATE LOCATION OF RED UNITS. 

C 

OISMAX-5000.0 
67 DO 9 I-l.NRU 

IF (lUSTAT (I) .EQ.2) GOTO 9 
IF (lUSTAT (I) .EQ.O) GOTO 76 
NF (I) -NF (I) '►I 
IF (NF (I) .LT.NOO) GOTO 9 
NF (I) -1 

76 DO 11 J • 1. NRU 

IF U .EQ. I) GO TO 11 
IF (lUSTAT (J) .EQ. 2) GO TO 1 1 
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DIST - X (1) - X (J) 

IF (DIST .GT. DISMflX) GO TO 9 
11 CONTINUE 

11 m - II m ♦ 1 
K7-II (I) 

X (I) -XIC (I,K7) 

Y m »Y1C (1.K7) 

MVTDIR(I)=»IDIR(I,K7) 

IPRDIR(I)-IDIR(I,K7) 

9 CONTINUE 
C 

C LINE— OF-SIGHT CHECK BETWEEN UNITS AND TARGETS SELECTION 

C 

DO 17 J-K,L 
NT U1 -0 

17 CONTINUE 

DO 12 1«1,NRU 
NT (I) »0 

IF (lUSTAT (1) .EQ.2) GOTO 12 
DO 16 J»K,L 

IF (lUSTAT (J) .EQ.2.0R.IUSTAT (J) .EQ.3) GO TO 16 
XX1«X (I) 

YTl-Td) 

133 CALL ELEV(XX1,TT1,TMAC1) 

XX2«X CJ) 

TY2-Y (J) 

CALL ELEV(XX2,TY2,TMACJ) 

LATQB-1 

LBTOA-1 

CALL LOS (XX1,TT1,TMACI ,0.0,SIZETK,XX2, YY2,TMACJ,0.0,S1ZETW. 
1LAT0B,LBT0A,V1SFRA,VISFRB) 

VISFR (1, J) -VISFRA 

VISFR (J, I) -VISFRB 

IF (VISFRA.GT.2L) GOTO 18 

LOST (I, J)«0 

LOST U, n »Q 

NLOSC (I , J) =NLOSC (I . J) 

NLOSC (J, I) «NLOSC (1, J) 

GOTO 16 

18 LOST (I, J)*l 
L0ST(J,I)-1 
NLOSC (1, J) »0 
NLOSC (J.D-O 

RANGE»SQRT ( (X (I) -X U) ) (Y (1) -Y (J) ) «m2) 

IF (RANGE. LT.RMINTK. OR. RANGE. GT.RMXTK) GOTO 20 

IF (Q(I, J) .EQ.1.0) GOTO 20 

lUSTAT (I) »1 

NT (I) =NT (I) +1 

M=NT (I) 

LOT (I,M) »J 
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ROT (1,M) -RANGE 
IF (H.EQ. 1) GOTO 20 
CALL SORT (1,M) 

20 IF (RANGE. LT.RMINTW. OR. RANGE. GT.RMXTW) GOTO 16 
IF (Q(J»n .EQ. 1 . 0 ) GOTO 16 
lUSTAT (J) -1 
NT (J) -NT (J) +1 
M-NT (J) 

LOT (J,M) -I 
ROT (J,M) -RANGE 
IF (M.EQ. 1) GOTO 16 
CALL SORT(J»M) 

16 CONTINUE 
12 CONTINUE 

00 25 1 = 1, NRU 

IF (lUSTAT (I) .EQ.2) GOTO 25 
IF (NT (I) .NE.O) GOTO 25 
lUSTAT (11-0 
NF (1) =0 
25 CONTINUE 
DO 79 J»K,L 

IF (lUSTAT (J) .EQ.2. OR. lUSTAT (J) .EQ. 3) GO TO 79 
IF (NT (J) .EQ.O) lUSTAT (J) =0 
79 CONTINUE 
C 

C UPDATE OF THE ACCUMULATED DETECTION PROBABILITIES. 
C 

lAA-1 

IBB-NRU 

ICC-K 

IOO«L 

FR-TOWFR 

OP»PM 

37 DO 14 I=IAA,1BB 

IF (lUSTAT (1) .EQ.2. OR. lUSTAT (1) .EQ. 3) GO TO 14 

DO 19 J-ICC.IDO 

PROP-0.0 

IF (lUSTAT (J) .EQ.2. OR. lUSTAT (J) .EQ. 3) GO TO 19 
OLDQd, J)=Q(1, J) 

IF (LOST (1, J) .EQ.O) GOTO 15 
IF (NT (I) .GT.O) GOTO 22 
PCTVlS=VISFR(J.n 

CALL LAMDAd, J.PCTVIS.DETRAT.PSUBK) 

QV=EXP (- (FL (I) >«OETRAThOP«DELThFL (J) ) ) 

IF (NT (J) .GT.O) GOTO 23 
Q d. J) =Q (I, J) mQV 
GOTO 19 

23 QP» (1 .0-PSUBK) «« (FRnOELTmFL (I) «FL (J) ) 

Qd. J) =Q(I, J) « (QV-^QP-QVnQP) 

GOTO 19 
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22 N5»NT (IJ 

00 24 I1»1,N5 
K1=L0T (Ml) 

flNGl»flTflN2 (Y (Kl) -T (I) ,X (Kl) -X (in 
flN62=flTflN2 (Y (J) -Y (1) ,X (J) -X (in 
IF ((fiNGl«flNG2) .GE.0.0) GOTO 77 
IF (flNG2.LT.0.0) GOTO 32 
flNG=2HPflnflNGl-flNG2 
GOTO 35 

32 flNG*2«PflnflNG2-flNGl 

35 IF (flNG.GT.Pfll) flNG=2xPflI-flNG 
GOTO 33 

77 flNG»flBS (flNG2-flNGl) 

33 fla=l5.0HpflI/l80.0 

IF (fiNG.GT.flfl) GOTO 24 
PROP=PROP^PTT (inN5) 

24 CONTINUE 

IF (PROP. EQ. 0.0) GOTO 34 

IF (NT (J) .GT.OJ GOTO 36 

CALL LflMDfl (K J.PCTVIS.DETRRT.PSUBK) 

DETRflT=DETRRTHRF 

QV = EXP (- (FL (I) KPROpKOETRflTKDELTKFL (J) ) ) 

Q (I » J) (1, J) mQV 
GOTO 19 

36 Q (I , J) »0.0 
GOTO 19 

34 IF (Iftfl.EQ. 1) GOTO 19 
Q (I , J) =1 .0 

GOTO 19 

15 IF (NLOSC (I. J) .LE.3) GOTO 19 
Q(KJ)=1.0 
19 CONTINUE 
14 CONTINUE 

IF (Iflfl.EQ.K) GOTO 38 
FR»TNKFR 

IBB=L 
ICC=1 
IDO=NRU 
0P=1.0 
GOTO 37 
C 

C FIRE ALLOCATION. 

C 

38 DO 28 1=1, L 
28 NA (I) =0 
DO 26 1=1, L 

IF (lUSTAT (I) .EQ.2.0R. lUSTAT (I) .EQ.3) GO TO 26 
IF (NT (I) .EQ.O) GOTO 26 
DO 27 J=l,3 
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flPOfl (I . J) =0.0 
27 CONTINUE 

IF INT 11) .EQ. 1) GOTO 78 
IF (NT (1) .EQ.2) GOTO 29 
NOT-3 

HM1=L0T (1.1) 

MM2»L0T (1 ,2) 

MM3*L0T (1.3) 

PftOB= (1 .0-0 (1 .HHl) ) hQ (1 .MM2) xQ (1 .MM3) 
flPOfl (1 . 1) =flPOfl (1. 1) +PTT (1 . 1) xPPtOB 
PR0B=Q(1.MM1) k(1.0-Q(1.MM2) )mQ(I,mm3) 

(TPOfl (1.2) =ftPOfl ( 1 . 2) *PTT (1 . 1 ) xPROB 
PROB»Q (l.MMl) xQ (1.MM2) x (1 . O-Q ( 1 . MM3) ) 
flPOfl (1 . 3) xflPOfl (1.3) »PTT (1 . II xPROB 
PROB= (1.0-0 (1 . MMll ) X (1.0-0 (1 .MM2) ) xO (1 .MM3) 
flPOfl (1. 1) »flPOfl (1. 1) +PTT (1,2) xPROB 
flPOfl (1 . 2) =flPOfl (1,2) -►PTT (2.2) xPROB 
PROB= (1 .0-0 (1 ,MM1) ) xQ (1 . MM2) x (1 .0-0 (1 .MM3) ) 
flPOfl (1,1) =flPOfl (1,1) ^ATT (1.2) xPROB 
flPOfl (1 .3) -flPOfl (1 .3) -»PTT (2,2) xPROB 
PROB-0 (I .MMl) X (1 .0-0 (1 .MM2) ) X (1.0-0 (I ,MM3) ) 
flPOfl (1,2) =flPOfl 1 1 , 2) *PTT ( 1 . 2) xPROB 
flPOfl (1,3) -flPOfl (1,3) "ATT (2.2) xPROB 
PROB- (1.0-0 (1, MMl) ) X (1 .0-0 (1 .MM2) ) x (1 . 0-0 (1 . MM3) ) 
flPOfl (1 . 1) -flPOfl (1. 1) +PTT (1.3) xPROB 
flPOfl (1.2) -flPOfl (1,2) +PTT (2,3) xPROB 
flPOfl (1,3) -flPOfl (1,3) ♦PTT (3,3) xPROB 
30 00 44 J-l.NOT 
KK-LOT (1. J) 

Nfl(KK)-Nfl(KK) + l 
IN-Nfl (KK) 

LOfl(KK,lN)-l 
POfllKK, IN) -flPOfl (l.J) 

44 CONTINUE 
GOTO 26 
29 NOT-2 

MMl-LOT (1,1) 

MM2=L0T(1,2) 

PROB- (1 .0-0 (l.MMl) ) xQ (1 ,mm 2) 
flPOfl (1. 1) -flPOfl (1. 1) -PTT (1. 1) xPROB 
PR0B=0( l.MMl) X (1.0-0 (I, MM2) ) 
flPOfl (1,2) -flPOfl (1,2) +PTT (1 , 1) xPROB 
PROB- (1.0-0 (l.MMl) ) X (1 .0-0 (1 ,MM2) ) 
flPOfl (1. 1) -flPOfl (1. 1) »PTT (1,2) xPROB 
flPOfl (1.2) -flPOfl ( 1 , 2) +PTT (2,2) xPROB 
GOTO 30 
78 NOT-1 

MM1=L0T(1,1) 

PROB-1 .0-0 (I .MMl) 

flPOfl (1 , 1) -flPOfl (1. 1) ♦PTT (1.1) xPROB 
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GOTO 30 
26 CONTINUE 
C 

C ATTRITION. 

C 

SUMRaO.O 
SUHB‘0.0 
DO 40 Nl.L 

IF (lUSTRT m .EQ.2.0R. lUSTAT (I) .EQ.3) GO TO 40 
H6^NR (I) 

3UM=0.0 

IF (M6.EQ.0) GOTO 47 
00 41 J»1.H6 
H7=L0R (I. J) 

IF (M7.LT.K) GOTO 42 

1TTPE=2 

GOTO 43 

42 ITTPE-1 

43 RRNGE»SQRT ( (X (I) -X (M7) ) mh2+ (T (I) -Y (M7) ) kx2) 

IF (ITRIT.EQ. 1) GO TO 131 

CALL STOCHUTTPE, RANGE, RJI) 

GO TO 5000 

131 CALL ETK (ITYPE, RANGE. T) 

AJ1=1.0/T 

5000 SUH-SUM+AJlxFL (M7) xPOAtl, J) kDELT 
41 CONTINUE 
47 OFL (I) =FL (I) 

FL m =FL 11) -SUM 
IF IFL m .GT.ZU GOTO 46 
FL (1) -0.0 
lUSTAT (1) =2 
46 IF II .LT.K) GOTO 60 
3UMB=SUMB+FL (1) 

TPOL 11) = IFO 11) -FL ID ) /FO ID 
GO TO 40 

60 SUHR-SUMR^FL ID 

TPOL ID = IFO ID -FL ID ) /FO ID 
40 CONTINUE 
C 

C PRINT AND CHECK FOR BATTLE DETERMINATION. 

C 

1TIME=ICh10 
DO 57 1=K,L 

IF IIUSTAT ID .E0.2) GO TO 57 
00 58 Jxl.NRU 

IF IIUSTAT IJ) .EQ.2) GO TO 58 
CHECK-X ID -X IJ) 

AVD=SQRT I IX ID -X IJ) ) xx2+ IT ID -T IJ) ) x«2) 

IF lAVD. LT.BREAK.OR.CHECK. lt. 50.) GO TO 250 
58 CONTINUE 
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C 

C 

c 



57 CONTINUE 
GO TO 99 



COMPLETE BLUE MOVE 

250 DO 251 I=K*L 

IF (IflLT.EQ. l.OR. IMOVE m .EQ. ITEM) GO TO 6000 
IF (lUSTflT (I) .EQ.O) lUSTfiT m »3 
IMOVE (I) *IMQVE (I) •►I 
IF (IMOVE (I) .LT. ITEM) GO TO 251 
X (I) =XR (1) 

Y(l)*Tfl(l) 

IF (lUSTflT (I) ,EQ.3) lUSTflT (I) =0 

251 CONTINUE 
99 WRITE(6,112) ITIME 

112 FORMAT (////IX, 'TIME = 14 , 1 X, 'SECONDS '//) 

WRITE (6, 113) 

113 FORMAT (IX, 'UN1T',5X, ' X ' , 8X. 'T ' , 5X, 'FORCE LEVEL ',2X, 'STATUS', 
12X, 'L0ST-PCT',2X, 'TARGETS') 

DO 59 1*1, L 
N6»NT (I) 

IF(N6.NE.0) GO TO 48 

WRITE (6,264) I , X (1) , T (I ) , FL (1 ) , I USTAT ( I ) , TPOL (1) 

264 FORMAT (3X , 1 1 , 3X , F7. 1 , 2X , F7. 1 , 6X , F3. 1 , 9X , 1 1 , 6X , F5. 3) 

GO TO 59 

48 WRITE (6, 114) I , X (I) , T (I ) , FL (I ) ,IUSTAT(1) ,TP0L(1) , 

1 (LOT (I, J) ,J*1,N6) 

1 14 FORMAT (3X, 1 1 ,3X,F7. 1,2X,F7. 1 ,6X,F3. 1 ,9X, 1 1 ,6X,F5.3,3X,3 (1 1, IX) ) 
59 CONTINUE 

C 

C CHECK FOR BATTLE DETERMINATION. 

C 

I0T»0 

DO 53 1-1, NRU 

IF (FL (I) .EQ.0.0) GOTO 53 

lOT-1 

53 CONTINUE 
IF(IOT.EQ.l) GOTO 54 
WRITE (6,117) 

117 FORMAT (IX, 'kh« RED FORCE IS ELIMINATED. END OF BATTLE.') 

GOTO 66 

54 I0T*0 

DO 55 1=K,L 

IF (FL(I) .EQ.0.0) GOTO 55 
I0T*1 

55 CONTINUE 
IF(IOT.EQ.l) GOTO 65 
WRITE (6,118) 

118 FORMAT (IX, '««« BLUE FORCE IS ELIMINATED. END OF BATTLE.') 

GOTO 66 
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551 6000 WRITE (6, 119) 

552 119 FORMAT (IX, '««« DISTANCE BETWEEN FORCES IS TOO CLOSE. END OF BATTLE 

553 1 •) 



55M 




GOTO 66 


555 


65 


JC=IC>1 


556 




GOTO 67 


557 


66 


STOP 


558 




END 



60 



nc^ 5 H ’ ,x’tri5 per : 






"iiNtaa ijtMttio 



Itli.f) Iflite ooo» 



fi tTOfl 

«i 

01(1 



tu 

S 28 

£ 2 t 

i<at 

•«« 

ftfS 



1 

z 

3 

« 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

2p 

21 

22 

29 

2M 

25 

26 

27 

28 

29 

90 

31 

32 

93 

34 

35 

36 

37 

38 

99 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 



SUfiftOUTINE L0S(Xfl.Tfl.TMACfl,TMlCfi,SIZ£fl,XB,tB.TMflCB.TMlCB.SIZE8. 
-LflTOB.LBTOR.VlSFftfl.VlSFRB) 

C 

CflKHON /HILLS/ XC (100) , TC (100) • PEAK (100) . SX (1001, ST (lOQl .AHO IIOOI 
COHHON /HILLS/ SCALE (lOO) . TUQAHO (100) , TWOSCL (1001 , BASE 
COMHON /HILLS/ NHlLLS 

COMMON /COVER/ CXC (ISO) , CTC (l50) . CPEAK (ISO) , CPXX (l50) , CPTT (150) 

. COMMON /COVER/ CPXT (150) , NCVELS 

COMMON /COUNTR/KH, KHW, KV. KN, KGRS, KELL, KINT 
COMMON /GRID/ LST (10, lO) , NHL (10, 10) , LI STH (450) .KHREP (100) ,KTR£f 
COMMON /GRID/ LSTC (10, 10) ,NC (10, 10) .LISTC (400) .KCREP (1501 
DIMENSION IGX (100) ,1GY (100) , lEL (100) ,CS1 (100) ,C32(i00) 

DATA NGR10/10/,GS1ZE/1000./ 

C SUBROUTINE TO COMPUTE FRACTION VISIBLE FOR OBSERVER TARGET PAIRS 
YISFRA-K 
VlSFRB-1. 

XBA«XS*XA 

TBA-TB-YA 

IF ( (XBA.EQ.Q.) .AND. (TBA.EQ.O.J ) RETURN 
IFISIZEA^TMICA.LE.O.) GO TO 5l0 
. IF (SIZEB^TMICB.LE.O.) GO TO 510 

IF (TMICA.LT.O.) VI SFRA- 1 . Q+TMI CA/SIZEA 

IF (TMICB.LT.O.) YlSFRB-1 . 0-^TMl CB/SIZEB 

ZA-TMACA TMICA ♦ SIZEA 

ZB-TMACB ♦ TMICB ♦ SIZES 

KTREP-KTREP^l 

ZBA-ZB-ZA 

XBASQ-XBAmXBA 

TBASQ-TBAmTBA 

XTBA-XBAmTBA 

TU0XBA-2.MXBA 

THOTBA-2,«Y0A 

C COMPUTE GRID SQUARES CROSSED BY A TO B LINE 
NGRSQ-0 

IF(XBA) 110.95.100 
95 XBA-0.1 

100 ISCX— 1 

X1NC>GS1ZE/XBA 
GO TO 120 
110 1SGX«1 

XINC— GSIZE/XBA 
120 IF(YBA) 140,125.130 
125 TBA-0.1 

130 ISGY--1 

TINC-GSIZE/TBfl 
GO TO 150 
140 lSGT-1 

TINC— 6SIZE/YBA 
150 IX-1^IFIX(XB/GSIZE) 

IF (IX.GT.NGRIO) 1X»NGR10 
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IT-1'^IFIX(YB/GSIZE) 

IF (IT.GT.NGRIDl IY«NGRID 
XNEXT«GSIZEh (FLOAT (IX) ♦O.Sn (ISGX-K)) 

YNEXT»GSIZEk (FLOAT (lY) +0.5 k (ISGY-1.)) 

XSTEP= (XB-XNEXT) /XBA 
YSTEP-(YB-YNEXT)/YBA 
160 NGRSQ=NGRSQ^1 
IGX (NGRSQ) *IX 
IGY (NGRSQ) -lY 

IF((XSTEP.GT.l.) .AND, (YSTEP.GT.l.)) GO TO 200 
IF (XSTEP-YSTEP) 170,1B0,190 
170 IX=»IX +I3GX 

XSTEP«XSTEP+XINC 
GO TO 160 
180 IX-IX+ISGX 

XSTEP»XSTEP-^XINC 
190 IY*IY+ISGY 

YSTEP-YSTEP-^YINC 
GO TO 160 

200 KGRS-KGRS+NGRSQ 

C GRID SQUARE LIST NOW COMPLETE IN IGX. IGY WITH NGRSQ ENTRIES 
C 

C NOW FIND WHICH COVER ELLIPSES TOUCH THE A TO B LINE, 

C CHECK ELEVATIONS AT SI AND S2 FOR EACH SUCH ELLIPSE 
NELS-0 
CHTMAX-0. 

IF (NCVELS.EO.O) GOTO 270 
00 260 Kel. NGRSQ 
IX^IGX (K) 
lY-IGY (K) 

N-NC (IX, lY) 

IF (N.EQ.O) GO TO 260 
LS*LSTC (IX, lY) 

LEND-LSi-N-l 
DO 250 L»LS,LEND 
KELL = KELL-^1 
IC-LISTC(L) 

IF (KCREP (IC) .EQ.KTREP) GO TO 250 

KCREP(IC)=KTREP 

RX*XA-CXC (1C) 

RY-YA-CYC(IC) 

PXX»CPXX(IC) 

PYY«CPYY(IC) 

PXY«CPXY (IC) 

AA«PXXmXBASO*^PYYkYBASQ^PXY«XYBA 

BB«PXXmTHOXBAkRX*^PYY*<TWOYBAkRY + PXYk (RXxYBA+RYxXBA) 

CC-PXXhRXkRX+PYYmRYkRY-^PXYnRXmRY-1.0 

ARG*BBkBB-4, OxAAxCC 

IF (ARG.LE.O.) GO TO 250 

SQ»SQRT (ARG) 
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SI— (BB^SQ) / (2.0«flfl) 

S2« (SO-BB) / (2.0mRR) 

IF(Sl.GE.l.) GO TO 250 
IF(S2.LE.O.) GO TO 250 
IFISI.LE.O.) GO TO 510 
IF(S2.GE.K) GO TO 5l0 
C NOW CHECK LOS RT SI RND S2 
KlNT-KlNT+1 
CPK-CPERK(IC) 

XS»Xfi-^S2«XBR 
YS-TR-^S2HYBfl 
CALL ELEV (XS.YS.HTS) 

HTS-HTS-^CPK 

ZS-ZR-^S2nZBR 

IF (LRTOB.EQ.O) GO TO 210 

CRLL KQVER (ZR. TMRCB, SIZES. ZB. S2, HTS. ZS. Vl SERB) 
IF (VISFRB.LE.O.) GO TO 510 
210 IF (LBTQfi.EQ.O) GO TO 220 
3-1. -S2 

CRLL KOVER (ZB. TMRCR, SIZER. ZR, S, HTS. ZS. VISFRR) 
IF (VISFRR.LE.O.) GO TO 510 
220 XS-XR^S1mXBR 
YS-Tfl*^SlMYBR 
CRLL ELEV (XS. YS.HTS) 

HTS-HTS^CPK 

ZS-ZR-^SInZBR 

IF (LRTOB.EQ.O) GO TO 230 

CRLL KOVER (ZR. TMRCB. SIZES. ZB, SI. HTS. ZS. VI SERB) 
IF (VISFRB.LE.O.) GO TO 510 
230 IF (LBTOR.EQ.O) GO TO 240 
S-l.O-Sl 

CRLL KOVER (ZB. TMRCR. SIZER. ZR.S. HTS. ZS, VISFRR) 
IF (VISFRR.LE.O.) GO TO 510 
240 NELS-NELS-^l 
lEL (NELS) -IC 
CSl (NELS) -SI 
CS2 (NELS) -S2 

IF (CPK.GT.CHTMRX) CHTMflX-CPK 
250 CONTINUE 
260 CONTINUE 
C RLL ELLIPSES CHECKED 
C 

C NOW STRRT ON THE HILLS 
270 DO 600 K-l.NGRSQ 
IX-IGX (K) 
lY-IGY (K) 

IF (NHL (IX. lY) .EQ.O) GO TO 600 
LS-LST (IX, lY) 

LEND-LS+NHL (IX. lY) -1 
DO 500 L-LS.LEND 
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I-LISTH(U 

IF(KHREPd) .EQ.KTREPJ GO TO 500 
KHREPm-KTREP 

C PROCESSING FOR HILL I STARTS HERE 
KH»KH^1 

C COMPUTE W -TOP OF THIS HILL ALONG 0-T LINE 
CX-XBA/SX (I) 

CY-YBA/SY m 

OX- (XA-XC (1) ) /SX (I) 

DY- (TA-YC (in /SY (I) 

FQ-TWOSCL (I) M (CXhOX+CYmDY+RHO (I) n (CXnQY'^CTmDXJ 1 
GQ-SCALE (I) H (CXkCX+CYhCY^TWORHO (I) *CXkCT) 

IF (GO .EQ.O.) GO TO 500 
W— FQ/ (2.mGQ) 

IF (ABS(W) .GT.5.) GO TO 500 
FSQ=FQ«FQ 

EQ-SCALE (I) M (OXmOX+DYkOY + TWORHO (DmOX^DY) 

POWER-EQ-FSQ/ (M.hGQ) 

IF (POWER .LT. -3.) GO TO 500 
HHW-PEAK (I) hEXP (POWER) 

KHW-KHW-^1 

IF(HHW.LE.BASE) GO TO 500 
ZH-ZA’^W><2BA 

IF ((W.LT.O.) .OR. (W.GT.l.n GO TO 300 
IF (HHW .GE. ZW) GO TO 5l0 
CVHTW-0. 

IF (NELS.EQ.O) GO TO 300 
DO 280 M-1,NELS 

IF ( (CSl (M) .GE.W) .OR. (CS2 (M) .LE.Wn GO TO 280 
IC-IEL (M) 

IF (CVHTW.lt. CPEAK (1C) ) CVHTW-CPEAK (IC) 

280 CONTINUE 

IF ( (HHW-^CVHTW) .GE.ZW) GO TO 510 
300 IF (HHW+CHTMAX.LT.AHINI (ZA-SIZEA.2B-SIZEB) ) GO TO 500 
C IF HE GET TO HERE THEN NEED TO FIND LOWEST SIGHT LINE OVER HILU 
C WEWTON ITERATION A TO B GIVING VlSFRB 
IF (LATOB.EQ.O) GO TO 400 
KV-KV+1 
V-W 

HHV-HHW 

NCT-0 

FV»FQmV 

TW0GV-2.hGQ«V 

330 FCNV-ZA-^HHVh (TWOGVhV + FV- 1 . ) 

KN-KN+l 

FACTOR- (TW0GVmTW0GV*^2.« (GQ + TWOGVhFQ) •►FSO) 

OFCNV-HHVhV»<FACTOR 

IF (BBS (DFCNV) .LT.l.E-10) GO TO 350 

V-V-FCNV/DFCNV 

IF(ABS(V) .GT.5.) GO TO 400 
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FV-FQmV 

TU0GV«2.mGQmV 

POWER - EQ^FV+GQmVmV 

IF (POWER .LT. -3.) GO TO 400 

HHV-PERK (I) xEXP (POWER) 

DHHV-HHVx (FQ-^TWOGV) 

ELV-ZR-^DHHVhV 

IF (RBS(HHV-ELV) .LT.l.) GO TO 350 
NCT-NCT^l 

IF (NCT.LT.IO) GO TO 330 
350 IF((V.LT,0.) .OR. (V.GT.l,) ) GO TO 400 

CVHTV»0. 

IF(NELS.EQ.O) GO TO 390 
DO 380 M»1,NELS 

IF ( (CSl (M) .GE.V) .OR. (CS2 (M) .LE.Vn GO TO 380 
IC«IEL (M) 

IF (CVHTV.LT.CPEfiK (IC) ) CVHTV = CPEfiK ( I C) 

380 CONTINUE 
390 HTV»HHV+CVHTV 
ZV»Zfi-^VxZBA 

CALL KOVER (Zfl, TMfiCB, SIZES, ZB, V, HTV, ZV, VI3FRB) 
IF (VISFRB.LE.O.) GO TO 510 
C NEWTON ITERATION B TO A GIVING VISFRA 
400 IF (LBTOA.EQ.O) GO TO 500 
KV»KV-^1 
V«W 

VMI»V-1. 

HHV-HHW 

NCT-0 

FV-FQmV 

TW0GV»2.«GQ«V 

430 FCNV-ZB-^HHVx ( (FQ^TWOGV) «VM1-1 .) 

KN*KN-^1 

FACTOR* (TW0GV«TW0GV + 2.m (GQ-^TWQGVnFQ) -^FSQ) 

DFCNV-HHVhVMImFACTOR 

IF (AB3 (DFCNV) .LT.l.E-10) GO TO 450 

V»V-FCNV/DFCNV 

IF(AB3(V) .GT.S.IGO TO 500 

VMl-V-1. 

FV*FQkV 

TW0GV«2.hGQhV 

POWER » EQ+Fv^GQhV«V 

IF (POWER .LT. -3.) GO TO 500 

HHV*PEAK(I) kEXP (POWER) 

DHHV*HHV« (FQ-^TWOGV) 

ELV*ZB+DHHV«VMI 

IF (ABS (HHV-ELV) .LT.l.) GO TO 450 
NCT»NCT*^I 

IF (NCT.LT.IO) GO TO 430 
450 IF ((V.LT.O.) .OR. (V.GT. l.n GO TO 500 
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CVHTV*0. 

IF (NELS.EQ.O) GO TO tl90 
DO 1180 M»1,N£LS 

IF ( (CSl (H) .GE. V) .OR. (CS2 (M) . LE . VH GO TO MSO 
IC«IEL (Ml 

IF (CVHTV.lt. CPE flK (ID ) CVHTV»CPEflK (IC) 
yao CONTINUE 
1190 HTV=HHV+CVHTV 
ZV»Zfl+VxZBfl 
S»-VM1 

CALL KOVER(ZB»TMflCfl,SIZEfl,Zfl,S.HTV.ZV, VISFRfl) 
IF (ViSFRfl.LE.O.) GO TO 510 
500 CONTINUE 
600 CONTINUE 
RETURN 

510 VISFRfl»0. 

VISFRB»0. 

RETURN 

END 
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SUBROUTINE SETUP 



C 

COMMON /HILLS/ XC (100) , YC (100) .PEAK (100) ,SX (100) ,ST (100) , RHO (100) 
COMMON /HILLS/ SCALE (1 00) , TWORHO (1 00) , TWOSCL (100) , BASE 
COMMON /HILLS/ NHILLS 

COMMON /COVER/ CXC (1 50) , CYC ( 1 50) , CPEAK (1 50) , CPXX (150) • CPTY (1 50) 

COMMON /COVER/ CPXT (1 50) , NCVELS 

COMMON /COUNTR/KH, KHW, KV, KN, KGRS , KELL , K1 NT 

COMMON /GRID/ LST (10, 10) , NHL (10, 10) ,LISTH (450) , KHREP (100) , KTREP 
COMMON /GRID/ LSTC(10,10) ,NC(10.10) ,L1STC(400) ,KCREP(150) 

READ (8,7) NHILLS 
READ (8, 47) BASE 
47 FORMAT (FIO. 4) 

7 FORMAT (16) 

17 F0RMAT(5F8.3,F6.4) 

DO 50 I»l, NHILLS 

READ (8,17) XC (I) ,YC (I) ,PEAK (1) , SX (1) , SY (I) , RHO (1) 

50 CONTINUE 

READ (8,37)LST 
READ (8, 37) NHL 
READ (8,7)NHT0T 

READ (8, 37) (LI STH (I) , I»1 , NHTOT) 

37 FORMAT (1615) 

READ (8,7)NCYELS 
IF (NCVELS.EQ.O) GO TO 65 
DO 60 I-1,NCVELS 

READ (8,27) CXC (I) ,CYC (I) ,CPEAK (I) ,CPXX (I) ,CPYY (I) ,CPXY (I) 

27 FORMAT (3F10.4,3E13.7) 

KCREP (I) --2147483600 
60 CONTINUE 

READ (8,37)LSTC 
READ (8,37) NC 
READ (8,7)NCT0T 

READ (8, 37) (LISTC (I) , I»1 , NCTOT) 

65 DO 100 I-l, NHILLS 
SX (I)-SX (I) nK625 
ST(I)-SY(I) nI.625 
XC(I)«(XC(1)-500.)«100. 

TC (I) » (YC (I) -930.) NlOO. 

TWORHO (I) =2.«RH0 (I) 

SCALE (I) — 1./ (2.H (l.-RHO (I) ««2) ) 

THOSCL (I) -2.HSCALE (I) 

KHREP(I) =-2147483600 

C ALL VALUES NOW IN METERS ON 0 — 10,000 GRID 
100 CONTINUE 

KTREP»-2147483600 

KH>0 

KHW^O 

KY-O 

KN-0 
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53 


KINT-0 


54 


RETURN 


55 


END 
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SUBROUTINE ROUTE 

COMMON /GRP3/ NBU, NRU. FL (6) , FO (6) , NOI (3) , XI C (3, 200) , TIC (3, 2Q0I , 
llOIR (3,200) ,flVSP, ISPD 

1, lUSTflT (6) , 11 (6) ,LOST (6.6) . VlSFRfl, VI SFRB, SI ZETK. 

1S12ETW.NT (6) .NF (6) . SRF . DI SMfiX. 

INLflSC (6.6) .VISFR (6.6) . RMlNTK. RMXTK. RMl NTH. RMXTH, OP. TQWFR, TNKFR» 
IPTT (3.3) .RF.PQfl (6.6) .flPOfl (6.6) .LOfl (6.6) .NR (6) .OFL (6) .POL (61 
DIMENSION XL0C(3.20) .rLQC(3.20) .N(3) 

IF (lSPD,EQ-y) DST»80.y63 
IF (ISPD. EQ. 3) DST = 67.053 
IF (ISPD. EQ. 2) DST = 53.6y3 
IF (ISPD.EQ.l) DST*40.232 
DO 300 1-l.NRU 
READ (5. 15) N(l) 

15 FORMAT (12) 

NL«N (1) 

DO 200 1N*2.NL 
READ (5.201) XLQCS.YLQCS 
201 FORMAT (F6.1. IX. F6.1) 

XLQC (1. IN) sXLOCS 
TL0C(1.IN)-TLQCS 
200 CONTINUE 

XLOCd.D-XlCd.l) 

TLQC (1. D-TIC (I. 1) 
lOlR(l.l) *0 
NL»N(1) 

NUM^2 

00 305 J>1.NL 

XL-XLOC (1. J^l) -XLQC (1. J) 

TL-TLOC(I.J+n-TLOC(I.J) 

DIST»SQRT (XLhh2+TL«h2) 

T-ABS (YL) 

Z-T/XL 
ANGL*ATAN (2) 

DEG>ANGLm 57.2958 
IF(J.EQ.l) GO TO 320 
XLN« (DST-EXTRA) hCQS (ANGL) 

D1ST» (DIST-^EXTRA) -DST 
YLN» (DST-EXTRR) hSIN (ANGL) 

XI C (I.NUM) »X1C (l.NUM-l) ^XLN^XLE 
1F(YL.GT.0.) GO TO 325 
YLN*-YLN 

325 YIC (I .NUM) aYlC (I . NUM- 1 ) ^YLN^ YLE 
1F(YL.GT.0.) GO TO 3M0 
1D1R(I.NUM)«-1F1X (DEG) 

GO TO 3U1 

3y0 IDIR(I.NUM) =IFIX (DEG) 

341 NUM«NUM-^1 
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320 XLN«D3TkC 03 (ANGU 
TLN»D3T*<3IN (flNGU 
IF ITL.GT.O,) GO TO 310 
TLN--YLN 

310 IF (DIST.LT.DST) GO TO 315 
XIC (I.NUM) -XIC (I,NUM-n -^XLN 
TIC (I,NUM)»YIC (I,NUM-n "^YLN 
IFCYL.GT.O.) GO TO 342 
IDIR (I,NUH) »-IFIX (DEG) 

GO TO 343 

342 IDIRd^NUH) *IF1X (DEG) 

343 DI3T=0I3T-D3T 
NUM-NUH+1 

GO TO 310 
315 EXTRfl=iOIST 

XLE»EXTRflxC03 (flNGL) 
YLE-EXTRAhSIN (flNGU 
IFCYL.GT.O.) GO TO 305 
YLE»-YLE 
305 CONTINUE 
300 CONTINUE 
RETURN 
END 
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SUBROUTINE LfiMOfl (1 . J, PCTV I S, OETRflT , PK)| 

C SUBROUTINE TO COMPUTE DETECTION RfiTE /(OETRAT) OF TARGET J 

C BT OBSERVER 1 GIVEN THE VISIBLE FRflCtlON (PCTVIS) . 

C 

COMMON /GRPl/ IPROIR (6) , ISECWO (6) pMVTOIR (6) ,X (6) , Y (6) ,SPD (6) 

TCFftCT-1.0 

ZER0L»0. 00001 

Pfll«22. 0/7.0 

7 D» (ISECWD (I) mPAI/i80.0) /2.0 

BBB* (1.0/ (2. On (SIN (0) -DnCOS (D)) )) 

IF (RBS(BBB) .LT.ZEROL) BBB-0.0 
RRR» (-BBB)nCOS(O) 

IF (RBS (RRR) .LT.ZEROL) RRR-0.0 

0TRNG-RTRN2 ( (T (J) -T (I)) , (X (J)-X(im 

P0«IPR0IR(IJnPRI/180.0 

IF ((PDnOTRNG) .GE.0.0) GOTO 1 

IF (PD. LT. 0.0) GOTO 9 

RNGLE»2nPRI^0TRNG-PD 

GOTO 10 

9 RNGLE»2NPR1*^P0-0TRNG 
10 IF (ANGLE. GT.PRI) ANGLE»2nPR1 -ANGLE 
GOTO 2 

1 RNGLE*R8S(P0-0TRNG) 

2 IF (ANGLE. GT.D) GOTO 3 
OUP-PD+D 
DLOW-PD-D 

ANGLFT«OTANG+ (15. OnPRI/1 80. ) 

IF (ANGLFT.GT.OUP) ANGLFT*OUP 
ANGLRT»OTANG- (1 5. nPRI /I 80. ) 

IF (ANGLRT.lt. OLOW) ANGLRT = DLOW 

PK-BBBn (SIN (ANGLFT) -SIN (AnGLRT) ) '►AARn (ANGLFT-ANGLRT) 

IF (PK.LT.0.0) GOTO 3 
IF(PK.GT.l.O) GOTO 5 
GOTO 8 

3 PK»0.0 
OETRAT^O.O 
GOTO 6 

5 PK»1.0 

8 RANGE»SQRT ((X (J) -X(I) ) n«2> (Y (J) -Y (I)) nn2) 
RR«0.001nRAnGE/PCTVIS 
T0ANG-ATAN2 ( (Y (I) -Y (J) ) , (X (I) -X (J) ) ) 

A0»MVTDIR(J)nPAI/180.0 
HORVEL»ABS (SPD (J) nSIN (TOANG-AO) ) 

DEN0M*l.y53*^TCFACTN (0. 5978^2. 1 88n (RRhh2) -0. 5038nH0RVEL) 

IF (DENOM. LE. ZEROL) OENOM»ZEROL 
DETRAT=«0. 003+ 1 . 088/0ENQM 
OETRRT=^OETRATnPK 

6 RETURN 
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SUBROUTINE ELEV (X , T. TMflC) 

C 

COMMON /HILLS/ XC (100) , YC (100) .PEAK (100) ,SX (100) ,ST (lOQ) ,RHQ (IQQ) 
COMMON /HILLS/ SCALE (1 00) , TWORHO (1 00) . TWOSCL (IQO) , BASE 
COMMON /HILLS/ NHILLS 

COMMON /GRID/ LSI (10,10) .NHL (10,10) ,LISTH(45Q) .KHREP(IOQ) .KTREP 
COMMON /GRID/ LSTC (10. 10) . NC (10, 10) , LISTC (400) , KCREP U5Q) 

DATA NGRID/IO/.GSIZE/IOOO./ 

C FUNCTION TO COMPUTE TERRAIN ELEVATION FOR GIVEN X, Y COORDINATES. 
ZMRX=BASE 
1X-1^1F1X(X/GSIZ£) 

IF (IX.GT.NGRID) IX=NGR1D 
lY-l + IFIX (Y/GSIZE) 

IF (lY.GT.NGRID) 1Y=NGRID 
IF (NHL (IX.IY) .EQ.O) GO TO 150 
LS=LST(IX,IY) 

LEND»LS-»NHL (IX, lY) -1 
DO 100 L>LS.LENO 
I«LISTH(L) 

QX- (X-XC (in /SX (1) 

QXSQ»QXmQX 

IF (QXSQ .GE.9.) GO TO 100 
QY- (Y-YC (in /SY (1) 

QYSQ*QY«QY 

IF (QYSQ .GE. 9.) GO TO 100 
QXY*TWQRHO (I) mQXnQY 
FACTOR*SCALE (I) (QXSQ-^QYSQ-^QXY) 

IF (FACTOR.lt. -3.) GO TO 100 
HT»PEAK(1)«EXP (FACTOR) 

IF (HT.LE.ZMAX) GO TO 100 
ZMAX » HT 
100 CONTINUE 
150 TMAC=ZMAX 
RETURN 
ENG 
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SUBROUTINE STQCH (I , RANGE , fl) 

C 

C SUBROUTINE TO COMPUTE STOCHASTIC ATTRITION COEFFICIENT 
C 

COMMON /GRP6/ ALPHA (6) 

COMMON /GRP3/ NBU. NRU, FL (61 . FO (6) , NOI (3J , XI C (3. 2001 , TIC (3,2QQJ » 
1I0IR(3,200J ,AVSP,ISPO 

I.IUSTAT (61 , II (61 .LOST (6.6J . VISFRA, VISFRB.SI2ETK. 

ISIZETW.NT (6) ,NF (6) .SRF.DISMAX, 

INLOSC (6.61 . VISFR (6.6) . RMI NTK. RMXTK, RMI NTN, RMXTW. OP. TOWFR. TNKFR, 
IPTT (3.3) .RF.POA (6.6) . APOA (6.6) .LOA (6.6) .NA (6) .OFL (6) , POL (61 
IF (I.EQ.2) GO TO 5003 
A-ALPHA (I) ( (1 ,0-RANGE/RMXTW) hm2) 

GO TO 5004 

5003 A-ALPHA (I) h ( (1 . 0-RANGE/RMXTK) hh 2) 

5004 RETURN 

END 
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C 

C 

c 



SUBROUTINE ETK ( I , RANGE, T) 



SUBROUTINE TO COHPUTE EXPECTED TIME TO KILL A TARGET. 

COMMON /GRP2/ TA (2) , T 1 (2) . TH (2) , TM (2) , TF I (2) » TF2 (2) , TF3 (2) , 
IP (2,6) ,PHH (2,6) ,PHM (2,6) ,PKH (2,6) ,TF (2) 

IF(I.EQ.2) GOTO 5 
TF (1)«TF1 (I) 

GOTO 6 

5 IF (RANGE. GT. 1000.0] GOTO 7 

TF (I) *TF1 (I) - (TFl (I) (lOOO.O-RANGE) / 1 000.0) 

GOTO 6 

7 IF (RANGE. GT. 2000.0) GOTO 8 

TF (1)-TF2(I)- ( (TF2(1)-TF1 (1) ) « (2000.0-RANGE) /lOOO. 0) 

GOTO 6 

8 TF (I) -TF3 (I) - ( (TF3 (I) -TF2 (1) ) m (3000. O-RRNGE) /IQOO.O) 

6 J» (RANGE^250.0) /500.0 
1F(J.GT.6) J=6 

T»TA(I)+T 1 (n-TH(l)^((TH(l)>TF(l))/PKH(I,J))-^((TM(I)+TF(in/ 
IPHM (1 , J) ) X ( (1 . 0-PHH (1 , J) ) /PKH (1 , J) '►PHH (I , J) -P (I , J)) 

RETURN 

END 
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SUBROUTINE SORT (I, Ml 

COMMON /GRP5/ LOT 16,6) , ROT 16,6) 

DO 19 J»1,H 

IF (ROT (1,M) .GE.ROT U, J) ) GOTO 19 
21 R-ROT (I, J) 

NN=LOT (1, J) 

ROT (I, J) *ROT (1,MI 
LOTn,J)=LOT(l,M) 

ROT (I,M) *R 
LOT (I,M) -NN 
19 CONTINUE 
RETURN 
END 
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SUBROUTINE KOVER (ZO, iMflCT. SIZET, ZT, S, HTS. ZS, VISFRT) 
C 

IFIS.EQ.O.) GO TO 2000 
IF (HTS.GE.ZS) GO TO 2050 
HEXT=ZO-^(HTS-ZO)/S 
EVIST-flMftXl (HEXT,THflCT) 

IF (EVIST.GE.ZT) GO TO 2050 
IF (EVIST.LE.ZT-SIZET) RETURN 
VIS= (ZT-EVIST) /SIZET 
IF(VIS.LT.VISFRT) VISFRT»VIS 
RETURN 

2000 IF (HTS.LT.ZO) RETURN 
2050 VISFRT*0,0 
RETURN 
END 
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APPENDIX C 



Definition of Variables in Com-puter Program 



ALPHA(I) 


= Initial attrition-rate coefficient for 
stochastic attrition module. 


APOA(I,J) 


= The average proportion of the 3 attacker 

of unit i allocated to fire on unit i. 


AVSP 


= Average speed of moving attacking units. 


EPSAK 


= Breakpoint distance between attackers and 
defenders. 


DISMAX 


= Maximum distance allowed between attacking 
units before the leading unit is delayed. 


DIET 


= The straight-line distance between two 
movement nodes inputed by the user. 


DST 


= The distance to be moved each time step 
by attacking units. - 


FL(I) 


= Force level of unit i. 


FO(I) 


= Initial force level of unit i. 


I ALT 


= Denotes whether the user desires alternate 
defensive positions or not. 


IC 


= Counts number of time units a defender 
has been moving. 


IDIR(I,J) 


= Direction of interval in i^^ route. 


IKI) 


- Interval index for unit i. 


IHOVE 


= Number of time units a defender is allowed 
for moving to an alternate position. 


IPRDIR(I) 


= Primary direction of movement for unit i. 


IRTE 


= Denotes whether user v/ants to input routes 
or not. 


IS 


= Random number seed used for stochastic attrition 


ISECWD(I) 


- Width of search sector for unit i. 
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ISPD 


= Input variable to denote user's desired 
speed for attackers movement. 


ITEM 


= Input variable denoting number of time 
steps allowed for defender's move. 


ITIME 


= Current time, in seconds, of battle. 


ITRIT 


= Input variable denoting whether attrition 
v/ill be stochastic or deterministic. 


lUSTAT(I) 


= Current status of unit i. 


LCA(I, J) 


= The number of the attacker of unit i. 


LOST(I,J) 


= Denotes \vhether line-of-sight exists 
between unit i and 3 or not. 


LCT(I,J) 


= The number of the target of unit i. 


MVTDIR(I) 


= Movement direction of unit i. 


N(I) 


= Number of nodes input ed by user for 
route i. 


NA(I) 


= Number of attackers of unit i. 


NBU 


= Number of blue units. 


riF(i) 


= Number of time units unit i is allowed to 
fire at the same location. 


NLOSCd, J) 


= Num.ber of continous time steps that 

line-of-sight does not exist between unit 
i and unit 3. 


NCI (I) 


= Number of intervals in the i^^ route. 


MRU 


= Number of Red Units. 


NT(I) 


= Number of targets of unit i. 


OFL(I) 


= Force level of unit i during previous 
time step. 


P(I,J) 


s t 

= Probability of 1 round hit by unit i 

in range band 3. 


PHH(I,J) 


= Probability of a hit following a hit 
by unit i in range band 3. 
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PHM(I, J) 


= Probability of a hit following a miss 
by unit i in range band j. 




= Probability of a kill given a hit by 
unit i in range band j. 


PI-I 


= The proportion of time a moving unit 
is searching for targets* 


POA(I, J) 


= The proportion of the attacker of 

unit i allocated to fire on unit i. 


POL(I) 


- Percent of unit i lost during the 
current time step. 


PTT(I, J) 


= Proportion of surviving firepower 
allocated to the target if there 

are 3 targets available. 


PANGS 


= Current minimum distance between attackers 
and defenders. 


PMINTK 


= Kinimura effective range for attacking 
v;eapon system. 


Ri-IINT’.7 


= Minimum effective range for defending 
weapon system. 


RMXTK 


= Maximum effective range for attacking 
weapon system. 


RMXT7/ 


= Maxim.um effective range for defending 
weapon system. 


ROT(I,J) 


a The range of the target of unit i. 


SIZETK 


= Size of attacking vehicle. 


SIZETV/ 


= Size of defending vehicle. 


TA(K) 


= Time to acquire a target for k^^ weapon 
system type (k = 1 , 2 ), 


TFl(K) 


= Time of flight to 1000m for k^^ v/eapon 
system type (k = 1 , 2 ), 


TF 2 (K) 


= Time of flight to 2000m for k^^ weapon 
system type (k = 1 , 2 ), 


TF3(K) 


= Time of flight to for k^^ wearon 

system type (k = 1 , 2 ), 
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TH(K) 


= Time to fire a round following a hit 
for weapon system type k (k = 1,2) ♦ 


TI(K) 


= Time to fire first round after target 
has been acquired for weapon system 
type k (k = 1,2). 


TM(K) 


= Time to fire a round follov/ing a miss 
for weapon system type k (k = 1,2). 


TNKFH 


= Firing rate for attacking weapon system. 


TOV/FR 


= Firing rate for defending weapon system. 


TPOL(I) 


= Total percentage of lost since battle 
began for unit i. 


VISFRCl, J) 


= The fraction of unit i seen by unit j. 


VISFRA 


= Fraction of unit A as seen by unit E. 


VISFRB 


= Fraction of unit B as seen by unit A. 


X(I),Y(I) 


= Coordinates of unit i. 


XA(I),YA(I) 


= Coordinates of alternate position for 
defender i. 


XIC(I,J) 
YIC(I, J) 


= Coordinates of the interval endpoint 

of the route for unit i. 


XL,YL 


= Distance added to previous interval 
endpoint for vehicle to move DST during 
a time step. 


XLOC(I,J) 
YLOCd, J) 


= Coordinates of the node inputed by 

the user for the route of unit i. 
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